function [C]=RenyiCost(mu,p,gamma_g,gamma_b,sigma)

C=-1*((1/(1-sigma))*p*log(gamma_g^sigma+(1-gamma_g)^sigma)+(1/(1-sigma))*(1-p)*log(gamma_b^sigma+(1-gamma_b)^sigma)-(1/(1-sigma))*log(mu^sigma+(1-mu)^sigma));
end

